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We investigate the predictions for fermion masses in the minimal realistic non-supersymmetric 
SU(5) model with the Standard Model matter content. The possibility to achieve b — r unification is 
studied taking into account all relevant effects. In addition, we show how to establish an upper bound 
on the ultraviolet cutoff A of the theory which is compatible with the Yukawa couplings at the grand 
unified scale and proton decay. We find A ~ 10 17 GeV, to be considered a conservative upper bound 
on the cutoff. We also provide up-to-date values of all the fermions masses at the electroweak scale. 



I. INTRODUCTION 

The hierarchy problem, unification of fundamental interactions, and fermion mass puzzle are some of 
the main motivations for physics beyond the Standard Model (SM). In particular, if we believe in unification 
of electroweak and strong interactions then the so-called grand unified theories (GUTs) represent the most 
natural extensions of the Standard Model. 

The first grand unified theory — the Georgi-Glashow SU (5) model [1] — was introduced in 1974. In that 
model each generation of the SM matter is unified in the 5 and 10 dimensional representations, and the 
minimal Higgs sector is composed of two representations: 5# and 24//. The Georgi-Glashow model is 
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arguably the simplest GUT. It is very predictive but it is certainly not realistic. Namely, one cannot unify 
the SM gauge couplings at the high scale, the neutrinos are massless, and a high-scale unification of Yukawa 
couplings of the down quarks and charged leptons contradicts experimental findings for the masses of those 
particles. 

Since the simplest SU(5) GUT is not realistic it requires appropriate modifications. There is a number 
of ways of doing that, but in order to preserve the predictivity of the theory those modifications should be 
minimal. With this in mind we have proposed in Ref. |Q] the simplest possible extension of the Georgi- 
Glashow model that is in agreement with experimental observations. Phenomenological and cosmological 
aspects of our proposal have been analyzed subsequently in Ref. Q30 . 

In this work we study in detail the Yukawa sector within the proposed framework, and define the upper 
bound on the ultraviolet (UV) cutoff A of the theory from the constraints imposed by proton decay lifetime 
measurements. We show that this UV cutoff depends on the absolute value of the tau lepton Yukawa 
coupling at the unification scale and not on the difference of the bottom quark to the tau lepton Yukawa 
couplings, as expected from bottom-tau unification. We find A < 10 17 GeV. Furthermore, as an essential 
ingredient in our study we perform an up-to-date analysis of all the fermion masses at the electroweak scale. 
These values are especially relevant for numerical studies of viability of various GUT models. 

The paper is organized as follows: In Section II we describe the minimal realistic extension of the 
Georgi-Glashow (GG) model. In Section III the predictions for Yukawa couplings at the unification scale 
are investigated. There we also list updated values of all the fermion masses at the electroweak scale. The 
UV cutoff of the theory is defined and evaluated in Section IV. Finally, we conclude in Section V. 

II. THE MINIMAL REALISTIC EXTENSION OF THE GEORGI-GLASHOW MODEL 

The Higgs sector of the minimal realistic non-supersymmetric SU(5) model [2] is composed of the fields 
24 H = (Eg, E 3 , £( 3) 2), £(3,2), S 24 ) = (8, 1,0) + (1, 3, 0) + (3, 2, -5/6) + (3, 2, 5/6) + (1, 1, 0), 15 H = 
($ a ,$ b ,$ c ) = (1,3, l) + (3, 2, l/6) + (6, 1,-2/3), and 5b = (H,T) = (1, 2, l/2) + (3, 1,-1/3), while 
the matter content remains the same as in the GG model. Here we use the SM (SU (3) x SU{2) x U(l)) 
decomposition to set our notation. The Lagrangian of our model includes all possible terms invariant un- 
der the SU(5) gauge symmetry, and accordingly includes higher-dimensional operators in order to write a 
consistent relation between fermion masses. (Influence of higher-dimensional operators on gauge coupling 
constants |4|,|5|] is assumed to be negligible.) The role of the 15h dimensional Higgs is twofold: it gener- 
ates neutrino masses through a type II see-saw mechanism [6], and contributes to the unification of gauge 
couplings. 
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The possibility to achieve unification in the present context has been investigated in Ref. 12j,|3I]. There 
we showed that there are only three fields which can help to achieve successful unification. Those are 
S3 C 24//, &a C 15//, and C 15//. We present in Fig. Q] the appropriate parameter space that 
generates successful unification of gauge couplings at one-loop. It corresponds to the region bounded by 




M GUT (KTGeV) 



FIG. 1: The whole parameter space is shown where we can achieve gauge coupling unification. The points PI, P2, 
and P3 define the allowed region, and the corresponding boundary values for the masses of <f> a , ${,, and £3 are shown. 



the lines of constant M$„ = 130 GeV, M<s> h = 242 GeV, and Af Sa = M z . 

In fact, the maximal value of Mqut at the two-loop level is somewhat larger than the one that corre- 
sponds to the benchmark point PI in Fig. Q] Namely, M GUT = 4.5 x 10 14 GeV for M Sa = M Sg = M z , 
M^ a = 1.1 x 10 4 GeV, Mg> b = 242 GeV and a^j T = 37.1. With this set of values we can establish an 
accurate upper bound on the proton decay lifetime. In a model independent way the proton lifetime t p is 
bounded by the inequality \M '• 

t p < 6 x 10 39 a^ 2 UT (My /10 16 GeV) 4 (0.003GeV 3 /a) 2 years, (1) 

where a is the matrix element, and My is a common mass of gauge bosons responsible for proton decay. 
For our purposes we set a = 0.015 GeV 3 [8] and identify My = Mqut- (The main source of uncertainty 
in Eq. (Q} comes from the matrix element a. For an up-to-date discussion on a see |9]. For a review 
on proton stability see [10].) Using the two-loop values of ogut and Mqut mentioned above we find 
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FIG. 2: One-loop gauge coupling unification at benchmark point PI in the Minimal Realistic SU(5) model (solid) in 
comparison with the SM (dashed). Notice the asymptotic free behavior above the GUT scale. 



^two-loops) <; i.4 x jq36 y ears [3|] Clearly, our model could be tested or ruled out at the next generation of 
proton decay experiments llli . 

Conversely, one can invert Eq. (O to establish a lower limit on the GUT scale; recall, Mqut = My. 
If we take r p (p — > ir°e + ) > 5.0 x 10 33 years [12] as experimental input, and assume o- G \j T ~ 37 the 
region below Mqut = 1.1 x 10 14 GeV is excluded by proton decay. This limit defines our benchmark 
point P3 in Fig.Q] Once we impose proton decay constraints on Mqut we also get an upper bound on the 
scalar leptoquark mass from the benchmark point P3: M$ 6 < 10 8 GeV. In addition, if we consider the most 
natural implementation of the type II see-saw mechanism (large M$ a ) the mass of the scalar leptoquark 
comes out in the phenomenologically interesting region O(10 2 — 10 3 ) GeV. 

We have analyzed the unification scenario at three benchmark points PI through P3 showed in Fig. [T] 
As we explained, these three points define the limits of the parameter space that yields unification of gauge 
couplings in the minimal realistic SU(5) model at one-loop. The particle content of the model also implies 
that the gauge coupling exhibits asymptotically free behavior between the GUT scale and the cutoff A. This 
feature is shown in Fig.|2]for the benchmark point PI. The value of cxgut at the scale M (A > M > Mqut) 
is given by: a G ]j T (M) = oi g \j T (Mgut) + m Mq UT • ^ n or( ^ er to ma ke any specific statements about 
the behavior of the gauge coupling at and above the cutoff the full structure of the more fundamental theory 
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that represents the ultraviolet completion of our model needs to be specified. As we comment towards the 
end, there exists well-defined underlying theory that reproduces our model below the cutoff A which also 
supports the asymptotic behavior of the SU(5) gauge coupling up to the Planck scale. 

We have also shown that the minimal realistic 577(5) predicts the existence of light fields. In the bench- 
mark point PI the Higgs field S3 and the scalar leptoquark &b are very light, while in the benchmark point 
P2 the lightest fields are <£ a and <!>{,. Finally in the benchmark point P3 there are two light Higgses <£ a and 
S3. We can conclude then that the minimal realistic non-supersymmetric SU(5) model could be poten- 
tially tested at the next generation of collider experiments, particularly at the Large Hadron Collider (LHC) 
at CERN. The possibility to explain the baryon-asymmetry in the Universe in this context has also been 



studied |3J]. See also Ref. H 1 311 for the possibility of probing second and third generation leptoquark param- 
eter space with the IceCube neutrino detection facility. Since in principle all those fields could modify the 
predictions for the Yukawa couplings at the unification scale we will study their effect in the next section. 

III. FERMION MASSES AND BOTTOM-TAU UNIFICATION 

We shall investigate predictions for fermion masses in the framework of the minimal realistic SU(5) 
model. The RGEs for the Yukawa couplings are summarized in Appendix B. In order to compute the 
values of the Yukawa couplings at the GUT scale we need accurate values for the fermion masses at the 
electroweak scale as initial conditions for the RGEs. As a first approximation it is sufficient to consider 
the third generation only, and neglect the Yukawa couplings of the first and second generations in the RGE 
evolution. For future applications however we shall update the values of the fermion masses at the Mz scale 
for all the three generations. 

As input value for the top quark mass we take the latest world average from the Tevatron Electroweak 



Working Group 11411 . and extract the top quark running mass assuming that this value corresponds to the 
pole mass: 

M t 4a a (M t ) f a s (M t )\ 2 „ 

+ - V ' +10.95 V ' +0{a\). (2) 



mt(Mt) 3 7T V 7r 

For the bottom quark mass we will adopt a conservative value at low energies m&(m;,) = 4.20 ± 0.10 GeV. 



This value is compatible through QCD evolution with the experimental measurement at higher energies 1 15], 
and agrees with most of the low energy determinations [12]. Then, to obtain the initial conditions for 
Eq. dB6l ) the top and the bottom quark masses are evolved to the Mz scale at three-loops |Q, [jjl with 



a s {M z ) = 0.1176 ± 0.0020 UJ]. We also extract the values of the other quarks from the PDG [12] and 



calculate them at Mz using three-loop RGE evolution in QCD accounting for the matching conditions over 



the heavy quark thresholds ll6Lll7ll . 

The physical lepton masses are taken also from the PDG ifl^j 
calculated at the Mz scale though the relation: 

a(M z 



The corresponding running masses are 



mi(M 2 



Mi 



1 



7T 



1 3, Ml 
1 + - In —^r 

4 Mf 



0(a 2 ) , 



(3) 



with a(Mz)~ l = 127.906 db 0.019. The input values for all the fermion masses and the corresponding 
values at the electroweak scale are summarized in Table H These values at Mz substantially defer and 
represent better reflection of our current knowledge of fermion masses from the values first evaluated in 
Ref. [18] and later updated in Ref. fl. 



TABLE I: Input parameters for the fermion masses and their values at the Mz scale. Capital M denotes pole masses, 
while m(/i) are running masses. 





input value 


running mass at Mz 


t 


M t = 171.4 ± 2.1 GeV 


m t {M z ) -- 


= 170.3 ±2.4 GeV 


b 


m b (m b ) = 4.20 ± 0.10 GeV 


m b (M z ) = 


= 2.89 ±0.11 GeV 


c 


m c {m c ) = 1.25 ±0.09 GeV 


m c (M z ) -- 


= 0.63 ± 0.08 GeV 


s 


m s (2GeV) = 95 ± 25 MeV 


m s {M z ) 


= 56 ± 16 MeV 


u 


m u (2GeV) = 2.3 ± 0.8 MeV 


m u {M z ) 


= 1.4 ±0.5 MeV 


d 


m d (2GeV) = 5.0 ± 2.0 MeV 


m d {M z ) 


= 3.0 ± 1.2 MeV 


T 


M T = 1776.99lg;26 MeV 


m T (M z ) = 


1746.45±g;|| MeV 


/' 


M^ = 105.6583692(94) MeV 


m„(M z ) = 


102.72899(44) MeV 


e 


M e = 0.510998918(44) MeV 


m e (M z ) = 


0.4866613(36) MeV 



Let us now study the predictions for fermion masses at the GUT scale. The leptoquark & b contributes 
at one-loop to the running of the Yukawa couplings for charged leptons and down quarks, while the field 
$ a modifies the RGEs for charged leptons. The field S3 does not couple to matter, only to the SM Higgs. 
It contributes to the renormalization of the mass and the couplings of the SM Higgs but does not modify 
the RGEs of the Yukawa couplings. The equations for the running of the Yukawa matrices Ye and Yd are 
given by: 

l 6vr 2 = y d (3% m + Y U Y}Y D 9(/i - M*J , 

16vr 2 ^ = Y E (/3l M ± \y}Y v 8(/i - M*J + \y$Y v 9(/z - , (4) 

where (3f M are the SM beta coefficients (see Eq. (|Blb ). The RGE of the up quark Yukawas is not modified 
with respect to the SM. The contributions of the fields <£ a and & b is due to the operator Yj,55 15# of 
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the Yukawa potential that generates the interactions Y,j9^ n lr, and Y^d^^L- ln the above equation the 
contribution of the field <£ a has been taken from Ref. [20]. Notice that usually in SU(5) theories there is no 
relationship between the Yukawa couplings of neutrinos and charged fermions. 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

Log 10 (^/1GeV) Log 10 (^/1GeV) 



FIG. 3: One-loop evolution of the masses of the third generation in the benchmark point PI of the minimal realistic 
SU(5) model (solid) and in the SM (dashed). 



Let us analyze the scenario where the contributions of the field <3? a and the leptoquark $5 have been 
neglected in Eq. (0]). Notice that at one-loop the mass of the field <3? a in all the parameter space is much 
below the natural value for the see-saw mechanism (M$ a = 10 13 — 10 14 GeV), and the Yukawa couplings 
for neutrinos are expected to be small. In Fig. [3] we show the one-loop evolution of the top, bottom and 
tau masses respectively in comparison with the SM for the benchmark point PI. Since in this case we have 
neglected the neutrino Yukawa contributions, the RGEs in the minimal realistic SU(5) model are the same 
as in the SM (Eq. (IB6I >). and the evolution of the charged fermion masses in that model is modified with 
respect to the SM only through the change in the evolution of the gauge couplings. Consequently, we obtain 
values of the bottom and tau Yukawa couplings at the GUT scale that are only 1-2% away from the SM 
prediction. A similar situation will happen for the other two benchmark points. Thus, independently of 
the benchmark point the Yukawa coupling of the bottom quark will lie at high energies below the Yukawa 
coupling of the tau lepton, and therefore it is not possible to achieve unification of Y T and Y\,. 

If we include the contributions of <£ a and/or the leptoquark <3?b in the running the difference between 
Y T and Y& at the GUT scale will be increased. The reason is that both contributions are positive, and the 
coefficient of the neutrino Yukawa coupling in the RGE for Ye is larger than that for Yd- This will result 
into a larger tau Yukawa coupling, and as we will see in the next Section into an smaller UV cutoff. In order 
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to obtain a conservative upper bound on the UV cutoff we will focus our analysis on the first scenario where 
the extra contributions of the fields <£ a and &b are neglected. 



IV. THE UV CUTOFF OF THE THEORY AND PROTON DECAY 



In this section we study the possibility to establish an upper bound on the UV cutoff of the minimal 



realistic 577(5) model. The relevant Yukawa potential up to order 1/A is defined by 



Wukawa = e ljklm Y ab Iff 5£ + 10j? 10* { ^f^5 n H + 10« Y%> 10^ 5 l H ^f^) 

+ 5i„ ioj? if 5 6i + 5^ lof k« 5 fcfc + 5 h 10jf' y f 5 &fe 



+ 5 ■ y (6) fk- 15 ij 4- 5 • Y (7) ^- g fT) 
i- o at i ab Obj io H -t- o ai j o6 °6J a ' 13 ' 



where i,j, k, I, m and n are the SU(5) indices, while a, b and c are the family indices. Once 24// gets a 
vacuum expectation value (VEV) (24//) = a diag(2, 2, 2, —3, —3) the GUT symmetry 577(5) is broken to 
the SM gauge symmetry. Then, the Yukawa couplings for charged fermions read 12111 : 

Y V = 4(Y + Y T ) - 12^(Y« + Y^ T ) - 2^{AY^ - Y^ T ) , (6) 
Y D = -Y^ + ^(3Y (4) - 2y( 5 )) , (7) 

Ye = -y(3) +3 ^(y( 4 ) + y( 5 )). (8) 

For neutrino masses, on the other hand, we find: 

r0\2 



Ml/ = Y^(5°)+Y^ { -^ 



(9) 



where (<5°) and (H°) are the VEVs of the neutral components of <£ a and H, respectively. The Yukawa 
couplings for charged fermions are diagonalized as follows: U T Yjj Uc = Y^ ia,& , D T Yd Dq = Y^ iag and 
Eq Ye E = Yf* 

Using the relevant relations for the Yukawa couplings we find: 

y B -y D = yffl 1 ^, (10) 
^Zttogut a 



where we use a = X/y/30 and My = Mqut = y| A gcur with g GUT = 4irotGUT- For the physical 
Yukawa couplings the relation reads: 

yf* = E%D*Y$**ri c E + ElY^E^t^-. (11) 
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If we require that the theory remains perturbative at the GUT scale: \ckk\ = \{EqY( 5 ' E)^/ ^/An < 1; 
hence the upper bound on A can be parametrized as: 



• , M GUT 
V a-GUT \Y T - a 6l Y^\ 

Here a k% = (E^,D*) kl (D' c E) ik . We have set k = 3 because in that case one finds the strongest upper 
bound. Notice that when a 3 * > (a 3 * < 0) and real, the maximal value of the denominator in Eq. (fT2l is 
Y T — Irf (Y T + Yfe). Since the GUT scale in our model is very low we have to be sure that it is possible to 
satisfy the experimental constraints on the proton decay lifetime. In particular, before defining the upper 
bound on the UV cutoff we have to determine which is the pattern of the a k% coefficients that is in agreement 
with the constraints set by nucleon decay. 

Let us discuss the relation between the upper bound on A and proton decay. It was pointed out in 
Ref. [7] that in order to find the upper bound on the proton decay lifetime in the context of the minimal 
realistic SU (5) model the following conditions have to be fulfilled: 

E C = DB 1 , D C = EB 2 , (13) 

where 

/ e ia J \ 



(14) 



e if3 i 
ye*» J 

In this case we have the maximal suppression for the proton decay channels. Now, using the above relations 
we get cu = (Bf D T Y^ E)u / y/47T and the upper bound on A (neglecting the phases) is given by: 



2 Mgut n „ 

A ^ \ x K? TFT • (15) 

OtGUT \Y T - Y d \ 

Notice that this is the bound which corresponds to the case o 3i > 0. By definition it is consistent with 
fermion masses and proton decay simultaneously. Therefore, using the values of the Yukawa couplings 
at the GUT scale, the most conservative upper bound on A is coming from the scenario when we have 
maximal suppression of the proton decay channels. It is important to say that usually the cutoff of a GUT 
model is evaluated from the difference between Y^ and Y T [22]. However, here we have shown that the most 
conservative upper bound on the UV cutoff, as given by Eq. ([TBI , is defined by the difference between Y T 
and Yd, and not by the departure from b — r unification. 

In the previous section we have already addressed the issue of numerical values of Yukawa couplings 
at the GUT scale. Since the value of the tau Yukawa coupling is almost independent of the benchmark 
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point, the actual bound on the cutoff simply depends on the ratio Mqut / y/otcuT- We hence summarize in 
TablellTlthe values of the masses of the fields $ a , and S3 at each benchmark point and the corresponding 
values of acuT, Mgut, Y t {Mqut), and the ultraviolet cutoff A. Notice that the most conservative upper 
bound on the cutoff of the theory at one-loop level is \ u PP er ~ 10 17 GeV, which could be identified with 



the string unification scale [23]. 



TABLE II: UV Cutoff for the three benchmark points. All the mass scales in GeV 



benchmark point 






Mx 3 


a GUT 


Mgut 


Y t {Mgut) 


A 


PI 


8.83 x 10 8 


242 


M z 


38.9 


2.53 x 10 14 


0.0098 


2.3 x 10 17 


P2 


130 


242 


1.69 x 10 10 


38.1 


1.11 x 10 14 


0.0098 


1.0 x 10 17 


P3 


6.92 x 10 4 


3.68 x 10 7 


M z 


38.7 


1.10 x 10 14 


0.0099 


9.8 x 10 16 



We briefly investigate how our analysis holds at the two-loop level. In Fig. [TJ we show our results for 
the scenario which corresponds to the benchmark point PI. To insure the proper inclusion of boundary 
conditions II24I1 at Mgut we set a^ 1 \ GJJT = a<^ T — Aj/(127r), where {A] , A2, A3I = {5,3,2}. (In 



addition to the central values given in TableUwe use CKM angles and phases from [ 12] and a^{Mz) = 
0.1176 ± 0.0020, a 2 (M z ) = 0.033816 ± 0.000027 and a x {M z ) = 0.016949 ± 0.000005 as our input.) 

Comparison between Figs. [2] and [4] shows that the GUT scale is slightly larger in the two-loop case. 
On the other hand, <£ a is significantly larger than in the one-loop case. And, we get Y T = 0.00986 and 
accordingly A = 3.1 x 10 17 GeV for the benchmark point PI, slightly above the one-loop result. 

Before we summarize our results, let us address the issue of the possible origin of higher-dimensional op- 
erators we invoke in Eq. ([5]). After all, they play decisive role in making our SU (5) model realistic. Namely, 
they affect the Yukawa sector in a way that simultaneously allows for realistic charged fermion masses and 
efficient suppression of the gauge mediated proton decay. In particular, they modify the GUT scale relation 
between Ye and Yd that basically rules out the GG model. And, they violate another prediction of the GG 
model, i.e., Yn = Yy, which would prevent one from completely suppressing the gauge boson mediated 
proton decay |2D. In fact, our considerations of these modifications and proton decay constraints resulted in 
the upper bound on the cutoff A. 

As we have shown, the cutoff comes out to be significantly below the Planck scale. This means that 
we cannot resort to the Planck scale effects to generate necessary higher-dimensional operators. It is then 
natural to ask for the credible renormalizable model that would effectively mimic the original proposal 
below the cutoff. To this end we observe that the most minimal renormalizable setup that yields the original 
model requires introduction of the following two matter pairs: (5,5) and (10,10). These pairs can clearly 
have gauge invariant mass terms above the GUT scale that can be identified with the scale A. Once these 



11 



fields are integrated out the effective model below A would have exactly the same features as the original 
model considered in this paper. For example, the relevant operators that eventually modify Ye = Yd 
relation are Ojl0j5^5 and 6j5j24ff5. (10,10) pair is needed to modify Y\j = Yy . Clearly, such a simple 
renormalizable realization of the model would also imply asymptotic freedom of the SU (5) gauge coupling 
between A and the Planck scale. 



V. SUMMARY 



We have investigated the predictions for fermion masses in the minimal realistic non-supersymmetric 
grand unified model with the SM matter content based on SU (5) gauge symmetry. We have shown that it 
is not possible to achieve b — r unification in this context since the extra contributions to the running of 
the Yukawa couplings are always positive and larger for Ye- We pointed out that the upper bound on the 
ultraviolet cutoff of the theory is A upper ~ 10 17 GeV which is consistent with the predictions for Yukawa 
couplings and the constraints coming from proton decay. In addition, we have provided up-to-date values 
of all the fermions masses at the electroweak scale. 
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APPENDIX A: RGES OF THE GAUGE COUPLINGS AT TWO-LOOPS 



The relevant equations for the running of the gauge couplings at the two-loop level take the form 

3 



dln/i ~ 2ir ai 



W + ^2E^ a 'W a jW + i a iW E Tr[a^]. (Al) 



32vr 3 



l=U,D,E 



The general formula for b, L and by coefficients is given in [25]. Besides the well-known SM coefficients we 
have: 

M M 

b?» 



^3 



w 







30 



5 



V / 



12 




FIG. 4: The gauge coupling unification at the two-loop level for central values of low-energy observables. The two- 
loop SM running is presented by dashed lines. Solid lines correspond to the benchmark scenario PI with £3, 
and $ a below the GUT scale. Vertical lines mark the relevant scales: Mj 3 = Mz, M<s> b = 242 GeV, M<p a = 
2.3 x 10 13 GeV, M GUT = 3.4 x 10 14 GeV and A = 3.1 x 10 17 GeV. 
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which we incorporate at the appropriate scales. The Cu coefficients are [26]: 



Cu 



( 17 1 3 \ 

10 2 2 

3 3 1 

2 2 2 

2 2 



V 



( 108 72 n \ 

25 5 u 

24 56 n 

5 3 U 





Obviously, Yukawa couplings enter the gauge coupling running at the two-loop level. Thus, one needs to 
run them as well at the one-loop level for consistency. We use the SM one-loop equations for the Yukawa 
couplings that can be found, for example, in Ref. [26]. 
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APPENDIX B: RGES OF THE YUKAWA COUPLINGS 



The beta functions for the evolution of the Yukawa couplings in the SM are given by 



SM 



0. 



SM 
D 



T-Gu + -{Y^Y v -YlY D ) 

t-g d + 1(yIy d -y^Yu) 



Pe 



T-G E + -Y ] E Y E 



(Bl) 



where 



T = TV (y^Ye + 3y£Yu + 3Y^Y D 



(B2) 



and 



(Gu\ 

Gd 
\G E J 



( il 2 8\ 

on 4 u 



20 4 

1 9 

4 4 

9 9 

4 4 



\alj 



It is convenient 112711 to define the following matrices: 



(B3) 



M E = Y%Y E , M D = YlY D , 



Mr 



(B4) 



whose diagonalization is given by 



E^M E E = di a g{Y 2 ,Y 2 ,Y 2 ) , 
D^MdDc = diag {Y 2 ,Y 2 ,Y 2 ) , 
ulMuUc = diag {Y 2 ,Y 2 ,Y 2 ) , 



(B5) 



The CKM matrix is defined through V CKM = uId c . Taking into account that AA^ = 1, for A = U,D, E, 



one can derive RGEs for the diagonal elements of these matrices: 



4tt 



dec 



u 



dm fi 2 



T-Gu + ^af 



-TW- 



CKM|2 D 
ij I a j 



4tt 



47T 



da? 
din fj, 2 

daf 
din fi 2 



D 



a? 



CKM 1 2 £/ 
y I 



(B6) 



where T = T/(4tt), G ; = G«/(4vr) and a[ = (y/) 2 /(47r), with i = 1, 2,3 the family index. 
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